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Abstract. We study the boundary-driven asymmetric simple exclusion 
process (ASEP) in a one-dimensional chain with long-range links. Shortcuts 
are added to a chain by connecting pL different pairs of sites selected randomly 
where L and p denote the chain length and the shortcut density, respectively. 
Particles flow into a chain at one boundary at rate a and out of a chain at the 
other boundary at rate /3, while they hop inside a chain via nearest-neighbor 
bonds and long-range shortcuts. Without shortcuts, the model reduces to the 
boundary-driven ASEP in a one-dimensional chain which displays the low density, 
high density, and maximal current phases. Shortcuts lead to a drastic change. 
Numerical simulation studies suggest that there emerge three phases; an empty 
phase with p = 0, a jammed phase with p = 1, and a shock phase with < p < 1 
where p is the mean particle density. The shock phase is characterized with a 
phase separation between an empty region and a jammed region with a localized 
shock between them. The mechanism for the shock formation and the non- 
equilibrium phase transition are explained by an analytic theory based on a mean- 
field approximation and an annealed approximation. 
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1. Introduction 

The asymmetric simple exclusion process (ASEP) has been widely studied in the 
past decades [1]. It is a nonequilibrium driven diffusive system of particles subject 
to the exclusion interaction. Despite its simplicity, the ASEP describes various 
nonequilibrium processes such as bio-polymerization [2], surface growth [3], traffic 
flow l4j[5], for example. Furthermore, the ASEP in one dimension is exactly solvable 
via the Bethe ansatz [5] and the matrix product ansatz j7] or direct solution of recursion 
relations [5]. The exact solution contributes to deeper understanding of fluctuation 
phenomena [HI I10| and nonequilibrium phase transitions [S] . 

Most studies on the ASEP have been performed on one-dimensional (ID) chains. 
On the other hand, there are many quasi-lD systems which involve long-range links. 
If a polymer chain folds randomly, there arise contacts between different polymer 
segments which arc far apart along a backbone [11]. Those contacts can be regarded 
as long-range links for a transport process on a network of polymer chains. A gene 
regulatory protein diffuses along the DNA chain to search for its target gene [T^]. It 
can make a long-range jump by dissociation from and reassociation with the DNA. 
Recent studies show that real world traffic networks have a complex structure with 
long-range links [131 [Tl]. In this respect, interests are growing in the study of the 
ASEP on complex networks [151 [TBI [T7] . 

In this work, we address the question: What is the transport capacity of a complex 
network for particles interacting via mutual exclusion? In order to contribute to this 
issue we consider the ASEP on a ID chain with long-range links, called shortcuts, for 
open boundary conditions. It will turn out that the shortcuts cause a drastic change 
in the phase diagram of the ASEP. 

We start with a brief review of the ID ASEP. The original model is defined on 
a ID lattice. The lattice sites are either occupied by at most one particle or empty. 
Multiple occupancy on a site is prohibited [exclusion interaction]. Particles may hop 
to the left and right with a bias to one direction [ASEP] . A closely related equilibrium 
process is the so-called symmetric simple exclusion process (SSEP) where particles hop 
in both directions with equal rates. For periodic boundary conditions, the model has a 
trivial steady state where every microscopic configuration is equally likely irrespective 
of the hopping bias [6]. However, the bias is relevant for the dynamic scaling behavior. 
In the context of growing interfaces, the SSEP belongs to the Edward- Wilkinson (EW) 
universality class [TS], and the ASEP to the Kardar-Parisi-Zhang (KPZ) universality 
class [H]. Both classes are characterized by the power-law scaling r ^ between 
characteristic time and length scales. The scaling exponent is given by z = 2 for the 
EW class and z = 3/2 for the KPZ class [Sll^. 

For open boundary conditions the chain is coupled to particle reservoirs at both 
ends. One acts as a particle source emitting particles at a rate a, and the other as 
a sink absorbing particles at a rate (3. Interestingly, the system driven by the open 
boundaries displays nonequilibrium phase transitions between low-density (LD), high- 
density (HD), and maximal-current (MC) phases [3 [5]. The system belongs to the 
LD phase if the capacity of the particle source a controls the particle fiux, i.e. if 
a < P and a < ph/2, where ph denotes the hopping rate of the particles. The system 
belongs to the HD phase when the outgoing rate /3 is smaller than a and pt/'^- The 
overall particle density is determined by the outgoing rate /3 and there is a macroscopic 
congestion of particles. When both a and /3 are larger than Ph/'^, the capacity of the 
system is limited by the capacity of the bulk. Such a phase is called the MC phase. 
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where the overall particle density is independent of a and /3. 

In the ID ASEP, particle hopping is a short-ranged process between neighboring 
sites. There are a few attempts to study the effect of a long-range hopping on the 
phase diagram of the ASEP for open boundary conditions. Szavits-Nossan and Uzelac 
considered the ASEP with a probabilistic long-range hopping [15] . A particle can hop 
to any site at a distance I with probability pi ^ l/l"'^^. They obtained that the phase 
diagram remains the same as that of the ID ASEP for <t > 1. For cr < 1, however, the 
system does not display any phase transition. Ha et al. considered a boundary-driven 
ID ASEP model where particles may perform a short-range hopping between nearest- 
neighbor sites or a long-range hopping [16]. Upon a long-range hopping, a particle 
jumps to an empty site directly behind a next particle in front of it. It was found 
that the long-range hopping introduces an instability towards a so-called empty-road 
phase. Otwinowski and Boettcher considered the ASEP on a one-dimensional chain 
decorated with hierarchically organized long-range links [TT]- This model displays 
the LD and HD phases. Besides, depending on the way the long-range hopping is 
implemented, an intermediate phase may also be realized. 

In this work, we investigate the role of long-range hoppings in a driven system 
in a generic setting. For that purpose, we study the ASEP on a ID chain with open 
boundary conditions and additional long range links connecting randomly-selected 
pairs of sites. The result is a graph similar to a small- world network [T3] . Additionally, 
we couple this network to two particle reservoirs, a particle source connected to one 
boundary site and sink to another. In contrast to the ID ASEP, particle source and 
sink now are connected via the small-world network. This setup describes a generic 
scenario for directed transport through a complex network with limited capacity. It 
is generally believed that long-range links make a system defined on a small-world 
network homogeneous [21]. In contrast, our study shows that the system develops a 
localized shock which separates a ID backbone into an empty region near the entrance 
and a fully-occupied region near the exit. Such an inhomogeneity is caused by the 
interplay between the boundary driving and the long-range hopping. The steady-state 
position of the shock depends on the particle input and output rates, which results in 
an interesting phase diagram. 

This paper is organized as follows: In Sec. [2l we introduce a boundary-driven 
ASEP in a ID chain with long-range links. In Sec. |3l we present numerical results 
obtained from Monte Carlo simulations. Measuring the overall particle density p, we 
obtain a numerical phase diagram which consists of an empty (E) phase, a jammed (J) 
phase, and a shock (S) phase. The overall particle density takes the value of p = (1) 
in the E (J) phase, while it varies continuously in the S phase. The shock phase is 
characterized by a localized shock which separates a ID backbone into empty and 
jammed domains. In order to understand the mechanism for the shock formation and 
the phase transition, we develop an approximate analytic theory. This is presented in 
Sec. m We summarize and conclude the paper in Sec. [SJ 

2. Boundary-driven ASEP with long-range links 

Consider a graph Q consisting of L sites which are labeled as i = 1, 2, • • • , X. Every 
pair {i,i + 1) with i = l,---,L— lis connected with a short-range link to form a ID 
backbone. In addition, we select pL pairs of sites at random and add long-range links 
between them (see Figure [Ij . The shortcut density, the total number of long-range 
links divided by L, is denoted by p. Every link has an orientation: A link between i 
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Figure 1. A graph consisting of L sites, short-range links, and long-range links. 
The site i = 1 and L are attached to a particle source and sink, respectively. 



and J is directed from i to j ii i < j and vice versa. Such a direction will be referred to 
as a forward direction. It is assumed that the boundary site « = 1, called an entrance, 
is attached to a particle source where particles are fed into the system at constant 
rate. The other boundary site i = L, called an exit, is attached to a particle sink 
which absorbs particles at a constant rate. 

We consider the ASEP process on such a graph. More precisely the particle 
dynamics is defined as follows: A particle enters the system at a rate denoted as 
a. Particles which have entered the system are then allowed to hop along forward 
links at rate one. If a particle is located at a site which is connected to others via 
multiple forward links one of these is randomly chosen. For each forward link the 
selection probability is given by 1/ki where ki denotes the number of forward links at 
a given site i. A particle can leave the system at site L (the exit) at a rate denoted 
as /?. Throughout the process, particles are subject to the exclusion principle which 
forbids multiple occupancy. Therefore any trial move violating the exclusion principle 
is rejected. 

One can represent the structure of a graph Q with an adjacency matrix A whose 
elements Aij take the value of 1 only if there is a forward link from site j to i. Due 
to the link directionality, Aij = for all i < j. Then, the hopping probability of a 
particle from site i to j is given by 



where = An denotes the number of forward links out of site i. It should be 
noticed that the hopping probability is a quenched random variable. It varies from one 
realization of a graph to another. Consequently, the quenched average over different 
realizations of the graph is necessary. 

Here, we are interested in the particle density distribution and the current. Let 
Ui (= 0,1) be the occupation number at site i. It is a stochastic variable evolving 
according to the ASEP dynamics. For a given realization of the graph, its average 
value is governed by the time evolution equation 




n. 



(2) 



for 1 < i < L and 




(3) 



J=2 




(4) 
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Here the angle bracket ((•)) represents the average of a quantity (•) for different 
reahzations of the stoehastic noise. 

In contrast to the ID-ASEP particle conservation in the bulk does not imply that 
the current is the same for every link. Therefore, the current can be defined in several 
ways. An obvious choice is to count the number of particles entering or leaving the 
system since the particle reservoirs are each connected via a single link. The incoming 
current Jm through the entrance and the outgoing current Jout from the exit are given 

by 

Jin = a(l - {ni)) , (5) 

Jout = fi{nL) ■ (6) 

In the bulk site i {— 1, ■ ■ ■ , L ~ I), the current Ji is defined as the total current of 
particles departing from sites j = 1, • • ■ , i and arriving at sites I — i + 1, ■ ■ ■ , L. It is 
given by 

i L 

Jr=Y^ ui-iin^il - ni)) . (7) 

With this definition of the current, the time evolution of the density is given by 
d{ni)/dt = Ji-i — Ji for 1 < i < L and d{ni) / dt = Jin — Ji aiidd{nL)/dt = J^^i — Jout- 
Hence, all values of the local currents J,;, Jm and Jout should be the same in the steady 
state, where d{ni)/dt = 0. 

The time-evolution equations ([2]), ([3]), and Q already illustrate the difficulty in 
calculating the quantities of interest. In order to determine the particle density two- 
point correlations have to be calculated, which in turn require higher order correlation 
functions. Furthermore, the quenched average over the random variables {u^} is 
necessary. An exact solution of this process is not available. So, we will investigate 
the model using a numerical simulation method in the following section. 

3. Simulation Results 

In this section we present numerical results obtained from Monte Carlo simulations. 
For a simulation, one generates a graph Q with an adjacency matrix A consisting of 
a ID chain of L sites and pL shortcuts. Particles on Q move in the following way: 
First, we select a random variable I G {0, 1, • • • , i} with equal probability, (i) If Z = 0, 
we add a particle to the entrance site i = 1 with probability a if it is empty, (ii) If 

< / < L, we try to move a particle (if any) at site I to a target site j selected among 
{I + I, ■ ■ ■ , L} with probability Uji given in ([1]). The particle move can be carried out 
if the target site is empty, (iii) li I = L, we remove a particle (if any) at the exit site 

1 = L with probability /3. The time is incremented by unity after (i + 1) trials. 

We are interested in the particle distribution and the current in the steady state. 
After a transient period of time interval Tj, we average the occupation number for the 
time interval to obtain the steady-state occupation number distribution (n^). Note 
that the mean occupation number varies from one realization of a graph to another. 
In a second step a quenched average [(• • ■)]g over graph realizations is necessary. The 
quenched averaged quantities will be denoted as 

P^ = [{n^)]g . (8) 
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Figure 2. (Color online) Particle density p (a) and current J (b) in the a/3 plane. 
The system size is L = 400 and the shortcut density p = 0.2. 



The overall particle density is given by 
1 ^ 



(9) 



The steady-state current is given by 

J = «(! - pi) = Ppl 



(10) 

The overall behavior of the particle density p and the current J is presented in 
Figure [5J Those data were averaged over Ns = 1000 graph realizations over a time 
interval Tg = 500000 after a transient interval Tt = 10000. Figure [2] suggests that 
there exist three different regimes. When a is small, the overall density p is close to 
zero. On the other hand, p is close to 1 when /3 is small. A finite-size-scaling (FSS) 
analysis reveals that the system indeed undergoes nonequilibrium phase transitions 
between three different phases, which will be discussed below. 

Figure |3lja) shows the plot of (1 — p) along a line a = 0.8 at several values of 
L = 100, ■ • • , 3200. Those data are weU fitted to a FSS form [22] 

(i-p)=i-"j:„((/3-/3e)i^/^") (11) 

with j3c — 0.15, ~ 0.4, and ~ 2.5. The scaling function has a limiting behavior 
^v{y ^ 1) ~ y^" with = XyV.^ ~ 1.0. The FSS indicates that system undergoes 
a continuous phase transition at a critical point /? = /3c- For /? < /3c, the density 
of vacant sites is zero. That is to say, the system is fully occupied by particles. 
Such a macroscopic state will be called a jammed phase. Near the critical point with 
€ ~ /3 — /3c , the density of empty sites scales as 

(l-p)^e^". (12) 

Figure [3Kb) shows the plot of p along a line /? = 0.8. The data are well fitted to 
a FSS form [22] 



with Ur 



P = 

0.36, Xp 



L-'''^Fj,{{a-a,)L^'''^) (13) 
~ 0.5, and Vp 2. The scaling function has a limiting behavior 



^p{y ^ 1) y p with /3p = XpVp ~ 1.0. Hence we conclude that the system undergoes 




Figure 3. (Color online) (a) Plot of the density of empty sites vs. jH at fixed 
a = 0.8. (b) Plot of the particle density vs. a at fixed /3 = 0.8. 




Figure 4. Phase diagram at p = 0.2. The location of symbols (□ between the 
E and S phases, and o between the S and J phases) is determined from the FSS 
analysis. The solid lines are the phase boundary obtained from the mean-field 
and annealed approximation developed in Sec. |3] 



a continuous phase transition at a critical point a = ac- The particle density vanishes 
for a < ac- Such a macroscopic state will be called an empty phase. Near the critical 
point with e = a — ac, the particle density scales as 

pr^e^'p . (14) 

Repeating the FSS analysis, we obtain the numerical phase diagram as shown in 
Figure SI The phase diagram consists of the empty (E) phase with p = and the 
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Figure 5. (Color online) Particle density profile in the shock phase. Each panel 
contains a data set obtained from a different realization of a graph. Monte 
Carlo simulation data are drawn with a red line. Also shown with a blue line 
are numerical data from the mean field approximation as explained in Sec. [J] 
Parameter values are L = 1000, p = 0.2, a = 1.0, and f) = 0.5. 

jammed (J) phase with p — 1. The other phase with < p < 1 will be called a 
shock (S) phase. Although the phase diagram looks similar to that of the ASEP on 
ID chains [3 [5], the nature of the phases is different. 

The particle distribution is intriguing in the S phase. Taking a = 1.0 and /? = 0.5, 
we have measured the particle density distribution {(rii)} in the steady state at a given 
realization of Q. Figure [S] shows typical distributions. There is a phase separation 
between a region with {ui) ~ and a region with {ni) ~ 1.0. The domain boundary 
between the two regions is called a shock. The localized shock gives a hint why there 
are phase transitions into the E and J phases. The system can be in the E (J) phase 
when the shock is absorbed at the exit (entrance). Hence it is crucial to understand 
the mechanism for the shock formation. This will be discussed in the following section. 

Before proceeding, we present a physical argument for the phase diagram. It is 
known that the diameter Id of the small-world network with nonzero shortcut density 
p scales as Id ^ ^ [SI- This means that any site can be reached from a given 

site within Id steps [131 [H]. Then, a particle injected at the entrance can move 
toward the exit with a diverging speed u~L/Z£)~L/lnL along the backbone in the 
L — > oo limit when a is so small that the exclusion is irrelevant near the entrance. 
Consequently, particles can escape instantly from the system and the E phase can be 
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realized. Using the same argument we can understand the origin of the J phase. A 
hopping of a particle from one site to another is equivalent to a hopping of a hole in the 
opposite direction. Hence, when (3 is small enough, a hole injected at the exit can travel 
toward the entrance with the diverging speed in OilnL) steps and the J phase can 
be realized. The S phase is a result of competition between the empty domain (which 
is stabilized by the long-range hopping of particles) and the jammed domain (which 
is stabilized by the long-range hopping of holes). This argument clearly shows that 
the small-world property is crucial in the formation of the E, J, and S phases. The 
small- world property emerges at any nonzero value of p |14] . Hence we expect that 
those phases replace the LD, HD, and MC phases of the conventional ASEP in ID 
chains immediately as one turns on the long-range links. 

We have also performed the numerical simulations with p ~ 0.1 and p = 1.0. In 
both cases we observe qualitatively the same phase diagram. As p increases, the E 
phase expands while the J phase shrinks. For example, when /3 = 0.8, we found that 
ac = 0.35, 0.36, and 0.39 at p = 0.1, 0.2, and 1.0, respectively. When a = 0.8, we 
found that /3c = 0.20, 0.15, and 0.07 at p = 0.1, 0.2, and 1.0, respectively. We do not 
find a simple explanation for the observed p dependence of the phase boundary. 



4. Analytic Results 

^.1. Mean field approximation 

We adopt a mean field approximation by assuming that 

{n,nj) = {n,){nj) (15) 
for all i and j. Applying the mean field approximation to ([2]), (|3]), and (|4]) yields that 

= + (l-(n,))Q, , (16) 

at 

where the auxiliary quantities are defined as 

\3>i 

and 

Q, = a5,s + u,jinj) | (1 - ^»,i) (18) 

\j<i 

with the Kronecker S symbol. The quantities Ri and Qi can be interpreted as particle 
evaporation and deposition rates at site i, respectively. They depend on the structure 
of an underlying graph through {uij} and the whole particle density distribution 
{{ni)}. This interpretation of the rates Ri and Qi allows to relate our model to 
the ID ASEP with constant evaporation and deposition rates which was studied in 
P51 [231 US . In contrast to that model the rates Ri and Qi in our model depend on 
the particle distribution as well as on the realization of the adjacency matrix. 
The steady-state density satisfying d{ni)/dt = is given by 

<"■> ^ fllf* ^ 

Since Ri and Qi depends on the particle distribution, (IT^ should be solved self- 
consistently. The self-consistent equation can be solved numerically via iteration. 
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Starting from any trial distribution, one updates it by evaluating the right-hand side of 
p^ . A particle distribution converges to a steady-state distribution without difficulty. 

In order to test the mean field approximation, we compare the particle density 
profile obtained from the Monte Carlo method and the mean field approximation. 
For a given graph Q, we have performed the Monte Carlo simulation and solved the 
self-consistent equation for the steady-state density profile. They are compared in 
Figure [5] The mean field approximation reproduces the phase separation and the 
shock in the density profile. In many cases, a mean field result is remarkably close 
to a Monte Carlo result (see the top and middle panels in Figure O . On the other 
hand, in some cases, there is a noticeable quantitative discrepancy between them (see 
the bottom panel in Figure [5)). Nevertheless, the mean field result still indicates the 
presence of the shock clearly. Therefore, we conclude that the mean field theory is 
suitable for the description of the phase transitions. 



4-2. Annealed network approximation 

The density profile given by the solution of ((T9)) depends on a graph realization 
Q. We have to perform the quenched average over graph realizations to obtain the 
disorder-averaged density profile pi = [{ni)]g. The quenched average is analytically 
intractable. Hence we further make an approximation by replacing an adjacency 
matrix element Aij, which is a quenched random variable, with its disorder- 
averaged value Such an approximation is called an annealed approximation. 
The annealed approximation is useful in studying physical systems on graphs or 
networks [26 l [27 l [28 | . 

On a graph of L sites, there are pL long-range links. Hence the probability 
to find a long-range links between two sites i and j 7^ i ± 1 is given by pi = 
2pL/{{L — l)(i — 2)) ~ '^p/L. Short-range links are connecting sites i and i + 1 
for i ~ 1, ■ ■ ■ , L — 1. Taking account of the short- and long-range links, we obtain that 

[^..]e='5,;,,+i + (l-<5„+i)-^ (20) 

for i > j and = for i < j. The parameter Uij is replaced by Ujj = [Aij]g / [kj]g 

with [kj]g = (1 + 2p{L - J - l)/L). 

In the annealed approximation, the self-consistent equation for the disorder- 
averaged density pi becomes 



Ri + Qi 

where Ri and Qi are given by 



(21) 




(23) 

with the boundary terms Rl = /3 and Qi = a. Equations pi|) . ([22|) . and ([23)) are the 
starting point for further analysis. 
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^.3. Shock state 

We first consider the sites with i = 0{1) near the entrance. Ignoring 0{L^^) 
corrections, one can approximate Ri and Qi as 

^" ^ nh;" -"•«>+ ilk" 

with the boundary term Q™ = a and the overall particle density p. These expressions 
allow for an interpretation for an effective dynamics near the entrance: A particle at 
site i performs a short-range jump to site i + l with probability 

W'J' = (26) 

or annihilates spontaneously with the probability 

W^r(p) = Y^(l-^) • (27) 

Due to the effective annihilation, the particle density should decay to zero 
exponentially with the distance from the entrance along the backbone with a 
characteristic length scale = 1/W^"(p) unless p = 1. This feature is consistent 
with the density profile shown in Figure [5j 

We next consider the sites i = L — I with I = 0(1) near the exit. Ignoring again 
0{L~^) corrections, one can approximate Ri and Qi by 

Br' = i-ft+i (28) 

with the boundary term i?^"* = P- These expressions suggest that particles near the 
exit have a following effective dynamics: A particle at site i performs a short-range 
hopping to site i + l with the probability 

wr' = 1 , (30) 

and particles are created spontaneously at each site with the probability given by the 
second term in (|29p . Due to the creation, the particle density should saturate to unity 
as one departs from the exit along the backbone. 

The density profiles stemming from the both boundaries converge to different 
values of and 1. So, there must emerge a shock as a domain boundary. The position 
of the shock along the backbone is denoted by is, which is related to the overall 
particle density as 

p=l-|. (31) 

Making use of the shock structure, the quantity Q""* in is given by 

nout _ , ^ 2p ^ 1 



L ^ 1 + 2p{L - ] - 1)/L 

3=is 



ln(l + 2pp) (32) 
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with 0{L ^) corrections being ignored. Correspondingly, the particle creation 
probability near the exit is given by 

W^r*(p)=Ml + 2w) . (33) 

So far we have established the shock state. Effectively, particles near the entrance 
hop to the right with the probability W™ or are annihilated with the probability 
W™{p). This dynamics results in a density profile {p™} which decays to zero as i 
increases. Near the exit, particles are created effectively with the probability W°'^*{p) 
and hop to the right with the probability W^"*. This results in a density profile 
which converges to unity as i decreases from L. Both profiles should be matched at a 
position ig = (1 — p)L to yield a shock. 

This situation is similar to the driven ID ASEP with particle creation and 
annihilation [531 [Ml Hi] or boundary driven multi-lane systems (see e.g. [211 150]). 
When the creation and annihilation rates are spatially uniform and inversely 
proportional to the lattice size, the system also develops a shock in the stead-state 
density profile [SSHMIES]. In comparison with the model studied in [23 ] [24 ] [25 ] [29 ] [30] . 
the creation and annihilation rates are not uniform in space: Particles are annihilated 
near the entrance and created near the exit. The difference results in the feature that 
the shock separates the empty and the fully jammed domains. 

The parameters W™{p) and W°'^*'{p) depend on the the overall particle density 
p = 1 — is/ L, which should be determined self-consistently. The overall density can 
be obtained from the current conservation, which will be explained in the following 
subsection. 



4- 4- Phase diagram 

Particles are injected at the entrance (i = 1), move to the right, and arc removed at 
the exit (i = L). So, the system can carry a nonzero current. The incoming current 
at the entrance and outgoing current at the exit are given by 

JM =«(i-pr) (34) 

Joutip) = PpT . (35) 

Because p™ and p^"* ^'I'C governed by the p-dependent effective dynamics, the incoming 
and the outgoing currents are given as a function of p. Particle number conservation 
requires that the incoming and outgoing currents should be the same in the steady 
state. The equality J™{p) ~ J°^^{p) determines the overall particle density p, hence 
the phase diagram. 

The effective dynamics is still too complex and does not allow for the closed- 
form solution for p™''°"*{p). Therefore, the exact phase diagram will be obtained from 
numerical solutions of the self-consistent equations and the current-balance condition. 
Before doing so, we apply an approximate scheme to the self-consistent equations in 
order to gain a physical insight. 

In terms of the effective dynamics, all particles introduced at the entrance have 
to be transferred toward the exit and therefore annihilated from the entry area. So, 
the incoming current can also be written as 

i>l 
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Using the explicit forms given in ([^^ and (|25p. the self-consistent equation for the 
density becomes as 



„2n 



'^^ 1 + 2p(l - p) + (p- 1 - 

for i>2. Because of the continuous annihilation of particles, we expect that p*" decays 
monotonically and rather fast to zero. In order to gain a qualitative understanding 
we can ignore — in the denominator to obtain that 

in 

^.n Pi _ (38) 



-(l + 2p(l-p)) 

Inserting these approximate solutions into p6p . one obtains that 

^ l + 2p(l -/>) ,„ , . 
1 + 2p ^1 ■ ^^^^ 



Comparing the two expressions for Ji„ given in ([34]) and p9p , one finds a solution for 
which yields that 

(1 + aj(l + 2p) - 2pp 

It is a decreasing function of p with Jm(0) = and J.m(l) = i+a(i+2p) • 

One can carry out a similar analysis to obtain an approximate expression for 
Jout{p)- First, the outgoing current given in ([55]) should be equal to the total particle 
creation rate, that is to say, 

jout = wr'Y.^i-p7') ■ (41) 

i<L 

It is easy to show that the void density (1 — p°"*) satisfies self-consistent equations 

1 _ out 

1 - 0°"* = f42) 

l + ln(l + 2pp) + (p°«* -p°fi) 

for i < L. Again the transport of particles via long ranged links can be understood as 
spontaneous creation of particles at sites close to the exit. Therefore, we expect that 
(1 — (0°"*) decays to zero as i decreases from L. So we can ignore (/0°"* — p°+i) in the 
denominator as in the previous case. A similar calculation then yields 

/3(l + ln(l + 2pp)) 
- l + /3 + ln(l + 2pp) ■ 

This is an increasing function of p with JoMt(O) = and JoMt(l) = i+/j+"n(i+2'pj • 

The particle number conservation requires that the incoming current and the 
outgoing current should be the same in the steady state. Figure [6] shows schematic 
plots of J™(yo) and J°"*(p) in three different situations. In the first case, the current 
curves may intersect with each other at p = po with < po < 1 as shown in Figuretn^a) . 
This case corresponds to the shock phase. The intersection determines the shock 
position ^5 = (1 — po)L. Ifp becomes greater (smaller) than po due a stochastic 
fluctuation, then the incoming current becomes smaller (greater) than the outgoing 
current. Consequently, the particle density is attracted toward the steady-state value 
and the shock is driven toward the steady-state position. This explains the reason 
why there is a localized sharp shock [331 [Ml [23 . 
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Figure 6. Schematic plots of Ji„ (solid lines) and Jout (dashed lines) against p. 

The second case with Ji„(p) > Jout {p) for all values of p is shown in Figure ID^b). 
Then, the shock is localized at the entrance and the steady-state density is equal to 
one. In this case the system belongs to the jammed phase. 

The third case is sketched in Figure [DJc) . If Jin (p) < J out (p) for all values of p, 
the shock is localized at the exit and the steady-state density is equal to zero. Then 
the system belongs to the empty phase. 

We have solved numerically exactly the self-consistent equations for p-" and p°"' 
to obtain the incoming and outgoing currents as a function of p at each set of values 
of a, (3, and p. The current balance condition allows us to evaluate the steady-state 
particle density and the current, hence the phase diagram. The resulting numerical 
phase diagram in the a/3 plane is presented in FigurelH The phase diagram a.t p = 0.2 
consists of the empty phase, jammed phase, and the shock phase, which is consistent 
with the Monte Carlo result. There is a quantitative discrepancy in the location of 
the phase boundaries, which is caused by the approximations. 

The annealed approximation requires the self-averaging property [311 132j . The 
self-averaging property in the small-world network was tested for the equilibrium 
Ising model |32j . To test the self-averaging property in our nonequilibrium model, 
we have measured the relative sample-to-samplc fluctuation of the particle density 

X = {[{I J2i r>'i)^h - [(r E» "»)]a)/[(r E» n^)]g■ in the E and J phases, it decays 
to zero as X ^ L^^^^ suggesting a strong self-averaging [31]. On the other hand, it 
converges to a finite value in the S phase, an indication of non-self-averaging. This 
tells us that the annealed approximation has a limitation. It explains successfully 
the mechanism of shock formation, but not its average position due to the strong 
fluctuations. A refined approach beyond the annealed approximation is necessary to 
study the sample-to-sample fluctuation phenomena in the S phase. 

5. Summary and Conclusion 

We have studied the boundary-driven ASEP in the ID chain with long-range links. 
This setup represents a generic situation for directed transport in a complex network, 
e.g., the exchange of data between two sites of a computer network. 

The backbone of nearest neighbor links ensures the existence of a path between 
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start and destination. The long ranged links add shortcuts to the transport network 
that are in principle able to enhance the capacity of the system. 

Considering the typical results for transport problems on complex networks one 
would also expect that long-range shortcuts, added randomly on to a lattice, are 
believed to suppress fluctuations and make a system homogeneous. This is, however, 
not the case when there is a boundary driving. Compared to the pure one-dimensional 
system it turns out that the long-range links play an essential role. They generate a 
localized shock which separates the ID chain backbone into an empty and jammed 
regions. 

Adopting the mean field and the annealed approximations, we have derived 
effective dynamics near the entrance and the exit, which are similar to those of the 
ASEP with spontaneous particle annihilation and creation, respectively. The effective 
theory reveals the mechanism for the shock formation and for the phase transition. The 
phase diagram consists of the empty, shock, and jammed phases. The shock phase is 
characterized by presence of a localized shock and separating a low and a high density 
domain. The shock position and the overall particle density p vary continuously with 
the model-parameters. In the empty (jammed) phase, the shock is anchored to the 
exit (entrance) to yield p = (1). 

In conclusion, our study shows that a driven system on a spatially disordered 
structure displays an inhomogeneous pattern. The appearance of a localized shock is 
reminiscent of one-dimensional systems without particle-conservation in the bulk. In 
contrast to these systems the mechanism driving the localization of the shock is neither 
the competition between bulk and boundary reservoirs nor an optimal partitioning 
between multiple lanes. Here, the long-ranged links enhance the mass transfer between 
entry and exit and thereby stabilize the position of the shock. 

Considering more generally the transport capacity of a complex network between 
arbitrary sites our results have important consequences. They indicate that, as far 
as the capacity of the feeding particle reservoir does not exceed the capacity of the 
exit reservoir, the sites of the backbone are only rarely occupied and can be used in 
parallel for transport issues between other sites. Contrary, overfeeding the backbone 
leads to a complete blockage of the sites in question which may spread over the whole 
network. 
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